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C " 3 ■ I review a recent analytic method for computing matrix elements of electroweak 

£f) ' operators based on the large-iV c expansion. In particular, I give a rather detailed 

description of the matching of the quark Lagrangian onto the meson chiral La- 
grangian governing K° — mixing. In the AS = 1 sector, I explain how to 
obtain an estimate for both e' /e and the AI = 1/2 rule. Finally, I give an example 
C D ' of how the method can also be useful for lattice calculations. 

■ , 1. Introduction 

Ph. 

The study of kaon weak interactions is difficult due to the fact that the 
strong interactions become nonperturbative at a scale Aqcd ~ 1 GeV 
which is higher than the kaon mass Mr-. However, because there is the 
. large hierarchy Aqcd/Mw <C 1 with respect to the scale of weak interac- 

| tions, Mwi the problem can be simplified with the use of Effective Field 

Theory techniques. In the Effective Field Theory of the Standard Model 
valid at the Kaon scale, one deletes from the Lagrangian any field whose 
associated particle mass is higher than Aqcd- At the same time one in- 
troduces higher dimension operators with the u, d, s quarks as degrees of 
freedom in order to reproduce the processes that in the Standard Model 
were taking place due to virtual exchange of the W and other heavy fields. 
The introduction of these higher dimension operators modifies drastically 
the ultraviolet properties of the Lagrangian making it look "nonrenormal- 
izable". However, if the basis of higher dimension operators is complete, it 
is possible to absorb all the divergences with the counterterms supplied by 
these operators. Of course, this takes care of the divergences, but it leaves 
a finite piece behind. Furthermore, it is a finite piece which is dependent 
on the conventions chosen to do the calculation (such as, e.g., the regular- 
ization scheme chosen -MS or MS-, the precise definition of 75 employed 
-NDR or HV-, evanescent operators, etc.). Obviously physical results 
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cannot be scheme/convention dependent and the ambiguity is resolved in 
the so-called matching conditions. These conditions equate a Green's func- 
tion computed in the Effective Field Theory and the full Standard Model, 
imposing that the physics is the same even though the heavy particles are 
missing in the former. 

The advantage of the Effective Field Theory technique is that it is sim- 
pler. This simplicity allows, e.g., the systematic resummation of all powers 
of the strong coupling a s accompanied by large logarithms in combina- 
tions such as a™ log™ Mw / Aqcd using renormalization group techniques 1 . 
However, below Aqcd there is no point in resumming powers of a s . Fully 
nonperturbative effects take place and, in fact, kaon transitions are de- 
scribed in terms of a Chiral Lagrangian with the kaon field as an explicit 
degree of freedom, and no longer in terms of fields for the u, d, s quarks. 
The matching condition between the quark Lagrangian and the Chiral La- 
grangian requires a nonperturbative treatment. This is where the large- iV c 
expansion comes in. This expansion is very well suited for this because it 
can be implemented both at the level of quarks and at the level of mesons 2 . 
However, even at the leading order in the large- N c expansion, Green's func- 
tions in QCD receive a contribution from an infinity of resonances whose 
masses and couplings are unknown. This is why one actually has to re- 
sort to an approximation to large- N c QCD. This approximation, which has 
been termed the "Hadronic Approximation" (HA) 3 , consists of the ratio of 
two polynomials whose coefficients are fixed by matching onto the first few 
terms in the chiral and operator product expansions of the Green's function 
one is interested in. The rational approximant so constructed constitutes 
an interpolator between the low- and high-momentum regimes, which one 
can use to perform the necessary calculations. 



2. AS — 2 operators and the K° — K° complex. 

2.1. Leading contribution: dimension-six quark operators. 

In the Standard Model, a double exchange of W bosons generates through 
the famous box diagram a AS = 2 transition amplitude between the K° 
and the K°. Below the charm mass, there arises the effective operator 

Q2 



c(u) l m X 2 c m 2 + m X 2 (m 2 ) + 2 % A c A t m 2 log^ }+h.c. , (1) 
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where ipL,R — ^2^^, \ — V* s Vid, m t . c are the corresponding quark 
masses, 771,2,3 are some numerical coefficients of 0(1) 23 , and 



c(» = (a s (»T 



1 + 



a 3 (n) /1433 1 
7T ^1936 + 8 



(2) 



expresses the running under the renormalization group of the operator in 
Eq. (1). The parameter n contains the scheme dependence and equals 
(—4) in the naive dimensional regularization (resp. 't Hooft-Veltman) 
schemes. In the case of the top one defines the effective mass 23 , 



(m 2 ) e// = 2.39M, 



w 



167 GeV 



1.52 



(3) 



to take into account that the top mass is heavier than the W. Physical 
observables such as the Kl — Ks mass difference and €k get a contribu- 
tion from the real and imaginary part, respectively, of the matrix element 
< K a \H s e f?\K° >. 

Below Kqcd ~ 1 GeV, it no longer makes any sense to think of an 
operator written in terms of quarks and, instead, one writes a Chiral La- 
grangian. The AS = 2 chiral operator at order 0(p 2 ) reads 5 



r S=2 _ 
L eff ~ 



G F 
16tt 2 



F 4 A 2 =2 Tr [A 32 (D^) UX 32 (D^) U] + h.c. , (4) 



where [A32]ij = 81382] is a (spurion) matrix in flavor space, F ~ 0.087 
GeV is the pion decay constant (in the chiral limit) and U is a 3 x 3 unitary 
matrix collecting the Goldstone boson degrees of freedom, which transforms 
as U — > RUL^ under a flavor rotation (R, L) of the group SU{2>) l x SU(Z)r. 



The scale A| =2 inherits some dependence on short-distance physics, 



■^S=2|dim-6 



9S=2 



2 2 2 / 2 \ 2 '^~^ J 't 

>7iA c m c + 772 \ {m t ) eff + 2-q s \ c \ t m c log — ^ 



iv- 



, (5) 



but there is also a coupling constant, gs=2, to be determined via a matching 
condition. 

What is this matching condition? Since covariant derivatives contain 
external fields I and r, i.e. D^U^U = d^WU + iWr^U - iP, one sees that 
Eq. (4) contains a "mass term", rf. s , for the external field . This is 

^ v ' 7 as M 7 as _ as 

precisely the external field that couples to the right-handed current d,R-f^SR 
in the kinetic term for the quark field in the QCD Lagrangian. Furthermore, 
a mass term for the r^ s field changes strangeness by two units, so that it 
can only come about from the quark Lagrangian because of the presence 
of the operator {sY^cLl) 2 in Eq. (1). Consistency demands that the two 
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mass terms be the same. Equating the term r^ s r^ s obtained from Eq. (4) 
to that obtained from Eq. (1) (plus gluon interactions) one obtains the 
matching condition 8 

9s=2 = c{n) 5s=2(a0 , where 

9S = 2 ^ ~ 1 ~ 32^if r(2-e/2) It W{zh (6) 

where n\ ad is an arbitrary scale used to define the integral in terms of a 
dimcnsionlcss variable z = Q 2 1 n\ a di where Q is a loop momentum. The 
coupling gs=2 is renormalization group invariant. The function W{z) is 
defined as 

W{z) = z^W^ RLR {z^ aA ), (7) 

where, in turn, W^ RLR is defined through an integral over the solid angle 
of the momentum q, 

J dn q g^Wfi&iq, OUactorizcd = " ff"") ^«L«(Q 2 ) ■ («) 

and Q 2 = - 9 2 . In Eq. (8) the function W^ R {q,l) stands for 

W^ v / R (q,l)= (9) 



,•3 ■ 



lim Lfixtfytfiz e^ +l y^(0\T{L^ d (x) R%(y) L" sd (0) <M}|0> 



i^o 
with 

L^ d (x) = s(x)r^ 1 d(x) , ^ s ( x ) = J(x) 7 ^i±^ s (x) . (10) 

Since W(z) ~ O(iV c ) and F 2 - 0(N C ) one sees that Eq. (6) yields 
9s=2 (m) = 1 — C(^c _1 )' with the unity stemming from the factorized part 
of the function W^ R " L ^ R . 

Notice that the function W^£ R {q, I) is essentially a 2-point "leftx left" 
Green's function, with incoming momentum Q, with a double insertion of a 
right-handed current at zero momentum (i.e. I — > 0). Even though W(z) is 
an order parameter of spontaneous chiral symmetry breaking and, therefore, 
receives no contribution from perturbation theory to all orders in a s , the 
integral in Eq. (6) is divergent (i.e. ill-defined) and must be regularized. 
Consistency demands that the regularization used be the same as in the 
calculation of the Wilson coefficient (2) and this is why the integral in Eq. 
(6) is done in MS in 4 — e dimensions. 
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In principle, knowledge of W(z) for the full range of z is required to 
calculate the coupling gs=2 (m)- However, lacking the solution of QCD at 
large N c , this information is not known. What is known, nevertheless, is 
the low- and high-z expansions of W(z) because they are given by chiral 
perturbation theory and the operator product expansion, respectively. If 
we build a good interpolator for the region in between, the answer is ready. 

In the large- N c limit the function W(z) is a meromorphic function, i.e. 
it has an infinity of isolated poles on the negative z axis, but no cut. The 
Hadronic Approximation (HA) I shall use is an approximation to this func- 
tion which consists in keeping only a finite number of poles, fixing their 
residues so that the coefficients of the chiral and operator product expan- 
sions are reproduced. In mathematics, this is called a rational approximant. 

In general the convergence of this type of approximants is difficult to 
establish for an arbitrary function 6 . However, in our case we expect it 
to work reasonably well. For one thing we expect W(z) to be a smooth 
function since it is a Green's function in the euclidean. Furthermore, if 
the operator product expansion sets in at scales > 1 GeV and the chiral 
expansion sets in at scales < M p , then the gap to interpolate by the approx- 
imant is not very large and, consequently, one may expect the error to be 
reasonably small. In particular, we have verified within a model that this 
approximation works nicely 7 and, in fact, the traditional success of vector 
meson dominance shows an indication that things may work similarly for 
QCD. At any rate, whenever possible, we shall check on the convergence of 
the approximation for the case at hand. 

Computing the chiral expansion for W(z) at low z one finds 11 ' 8 

W(z)t* 6-24 4^ (2L! + 5L 2 + L 3 + L 9 ) z +... , (11) 

where the first (resp. second) term comes from the contribution of 0{p 2 ) 
(resp. 0(p 4 )) in the strong chiral Lagrangian. 

At high z the operator product expansion yields 

(12) 

The leading term in 1/z in this expression stems from the contribution of 
dimension-six operators to the operator product expansion. The k param- 
eter is the same as that in Eq. (2) and contributes a nonvanishing result 
insofar as the calculation is done away from 4 dimensions, i.e. 4 — d = e =/= 0. 
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This is how our matching condition knows about short-distance renormal- 
ization scheme conventions. 
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Figure 1. Resonance diagrams contributing to W(z). 

The 5k term in Eq. (12) is the contribution from dimension-eight oper- 
ators in the operator product expansion. We may use this term to monitor 
the convergence of our expansion in the following way. First, we may start 
by considering the 5k term in the operator product expansion as "sublead- 
ing" and neglect its contribution altogether. After reinstating it later on, 
and redo the calculation, we may compare the two results thus obtained. 
Clearly, if our approximation is to make sense, both results should be close 
to one another. 

So let us neglect 5k in Eq. (12) for the time being. Since we wish 
to calculate the integral in Eq. (6), we need to interpolate between the 
low- and high-z regimes in Eqs. (11,12). To this end we observe that the 
pole structure of the diagrams depicted in Fig. 1 shows single, double and 
triple poles (those which can be obtained by cutting through the solid line 
between crosses). In the strict large- N c limit the Mittag-LefHer's theorem 9 
for meromorphic functions allows us to write W(z) in the form 

where the infinite sum extends over all possible resonances allowed by quan- 
tum numbers. The analytic structure of (13) suggests that, as a first step, 
we restrict ourselves to just one resonance in this sum. The HA interpolator 
so constructed reads then, 

(z + p) (z + p) 2 (z + pf 

where p = M 2 / ' n\ ad: and we determine the constants ay, by, cy by impos- 
ing that it reproduces the behavior in Eqs. (11, 12). The result of this 
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interpolation is shown as the dash-dotted curve in Fig. 2. This figure also 
shows the OPE and chiral expansion as dashed curves at large and small 
values of z, respectively. 




0.2 0.4 0.6 0.8 



z 

Figure 2. The dash-dotted curve corresponds to the interpolator W[z] in Eq. (14), i.e. 
without considering the Spc contribution, normalized to its value at the origin, versus 
z = Q 2 /fJ. 2 had for fi had = 1.4 GeV, F = 85.3 MeV and M v = 770 MeV. The solid curve 
is the full function W(z) in Eq. (19) after including the 5jf contribution in Eq. (12) (see 
text below). As one can see both curves are quite similar (the insert is a blow-up of the 
region z ~ 0.4). The two dashed curves correspond to the chiral and operator product 
expansions, for small and large values of z, respectively. 

Notice that the shape of the interpolating function W(z) shown in Fig. 
2 is very smooth. This kind of shape for the Green's function is very 
representative and we have found similar shapes in all the different cases 
we have studied. 

One can now use the example of W(z) to develop some intuition of the 
role played by the different energy scales in the problem. Although in a 
one-scale theory like QCD all scales are ultimately related to each other, 
it is a fact of life that there is a clear numerical hierarchy. First, there are 
chiral parameters such as or (ip'4>) 1 ^ 3 whose scale is ~ 100 — 200 MeV. 
Second, there is a typical resonance mass, or mass gap, whose scale is much 
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larger, i.e. Aqcd ~ 1 GeV. Notice that, since F v <~ Av 2 but Aqcd ~ A^? , 
a naive use of the strict large- 7V C limit would lead one to the erroneous 
conclusion that Aqcd/F^ is negligible. 

In matching conditions like (6) this observation is crucial. The reason 
is that the factorized contribution is governed by F n -which in Eq. (6) 
has been scaled out and this is why the first term is unity- whereas the 
unfactorized contribution is governed by Aqcd- 

Since W(z) ~ 1 in the region z < Aq CD / n\ ad -before the OPE sets in- 
and neglecting the logarithmic divergence from the OPE tail, this means 
that 



which says that the correction to unity in Eq. (6) will be of order 
Ag CD /(47ri^) 2 , which is not a negligible contribution at all. In this case 
keeping only the factorized contribution is not safe because the N c — > oo 
limit happens to select a scale which is "abnormally small" , F n , as compared 
to the larger scale Aqcd, which can only show up in the ncxt-to- leading 
1/N C terms. In these cases it is not unnatural to expect large unfactorized 
contributions. As a matter of fact, later on we shall see that the unfactor- 
ized contribution is ~ 50% of the factorized contribution in the case of Bk 
but it is several times larger in the case of the strong penguin contribution 
to e'/e. 

However, once one has included the unfactorized contribution with its 
Aqcd scale, there is no further scale in the game. Therefore, there is no 
reason to expect that subleading effects will still yield larger contributions 
with the consequent breakdown of the large-A^ expansion. These sub- 
leading effects are typically either an OZI-violating amplitude or related 
to resonance widths, and all evidence so far is compatible with these two 
being reasonably small. Notice that something similar to this happens in 
the relationship between the rf mass and the topological susceptibility for 
different number of colors 10 : even though the 77' mass is not at all small as 
compared to Aqcd -whereas it vanishes when N c — > 00- the l/N c expan- 
sion gives a good description of the lattice data for different values of N c , 
i.e. the expansion does not break down. 

Let us come back to our matching condition (6). Having disregarded 
the 5k contribution, one sees that the interpolator W(z) in Eq. (14) crosses 
the OPE curve at z ~ 0.4 (which is equivalent to Q ~ 900 MeV). Merging 
into the OPE only takes place at values of z which are larger than those 




2 



(15) 
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shown in the plot. Assuming that the OPE is a fair description starting 
from Q <~ 900 MeV onwards, one can now get an estimate of the integral 
in Eq. (6) by first integrating with the interpolator W(z) in the region 
< z < 0.4 and then with the OPE curve in the region 0.4 < z < oo. The 
result is usually presented in the form: 

B X K = \ gs=2 , (16) 

and one finds B\ — 0.38 ± 0.11, where the error is an estimate of higher 
order 1/N C corrections, ~ 30%. This error 8 covers all reasonable variations 
in the input parameters (e.g. the p mass). 

The combination appearing in Eq. (16) parametrizes the K° — K° 
matrix element of the AS = 2 operator in Eq. (1), 

< K°\ c(p) (s L ^d L f (p)\K° >= \f k M 2 k B k , (17) 

after including chiral corrections (this is why Bj( ^ B K ), and governs 
quantities such as Sk and the Kl — K$ mass difference. The inclusion of 
these chiral corrections is a necessary ingredient in order to make contact 
with the physical world since, moreover, there are indications that they 
may be sizeable 12 . We plan to be able to report on this in the near future 
but, for now, the result after Eq. (16) is still in the chiral limit. 

One may now include the contribution from the 8k term in the opera- 
tor product expansion of Eq. (12). This coefficient 8k parameterizes the 
following matrix element 

< 0\g s s L G« u \ a ^d L \K(p) >= -iV2F S 2 K Vli , (18) 

and a sum rule analysis 13 gives 8 K = 0.12 ± 0.07 GeV 2 . Since we now 
have one more condition at large z, we have to include one more resonance 
in the interpolator in order to achieve matching. First, we observe that 
the low-z behavior of W(z) depends on L 3 and this low-energy constant, 
unlike the other Lj's appearing in (11), receives contribution from scalars 
resonances 14 . Second, a simple pole contribution from a scalar resonance 
is compatible with the analytic structure and quantum numbers shown in 
Fig. 1. Therefore, we enlarge the interpolator in Eq. (14) with a scalar 
pole to read 

ttz/ \ a v W c v a s . . 

W(z) Vi S = - r - 7 + 7— ^- + 7— T3+7— 7 , 19 

(z + pv) {z + p v y (z + pv) 3 (z + ps) 
and allow for a generous variation of the scalar resonance mass Ms = 
900 ± 400 MeV. Imposing that W{z)y.s matches both expansions in Eqs. 
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(11, 12) one can obtain the unknown residues and determine W(z), Eq. 
(19). This function is the solid line plotted in Fig. 2. One sees that 
it matches the OPE in a smoother way than the interpolator obtained 
earlier (dash-dotted line in Fig. 2), but there is no dramatic difference in 
the area underneath. In summary, the neglect of higher order terms in the 
expansions (11, 12), as a first approximation, is a self-consistent procedure 81 . 

Using now this improved function W(z) and Eq. (6), we can calculate 
B\ once again and obtain 

B\ = 0.36 ± 0.15 , (20) 

as our final value. We emphasize that, because the coupling gs=2 in Eq. 
(5) is of 0(p 2 ) in the weak chiral Lagrangian (4), our value is in the chiral 
limit. This result is in nice agreement with recent determinations done 
on the lattice, where it now begins to be possible to consider dynamical 
fcrmions 15 . 

2.2. Corrections from dimension- eight quark operators. 

The passage from the Lagrangian in terms of quarks and gluons (1) to the 
Lagrangian (4) in terms of Goldstone mesons requires the integration of all 
the resonances. But in this integration the mass scale involved is the mass 
gap, i.e. Aqcd ~ 1 GeV, which is not negligible compared to the charm 
mass. Therefore, in the matching condition for A|_ 2 in Eq. (5) there may 
be extra contributions of 0(Aq CD ). 

To analyze the presence of these extra contributions one obviously must 
consider effects which are of 0(AQ CD /ml), relative to the contributions in 
Eq. (5). Therefore, after integrating out the charm quark, one must go 
to a Lagrangian of quark operators which are dimension eight rather than 
dimension six as in Eq. (1). Since in this case there is not the competition 
between the scales F n and Aqcd we discussed earlier, to simplify matters, 
we shall take the large- N c limit. 

Dimensional analysis shows that AS = 2 dimension-eight operators, 
unlike the dimension-six operator in Eq. (1), cannot come with a quark 
mass out front. It is not surprising, then, that the GIM mechanism becomes 
fully operational for scales m c < u < m t ,Mw while the charm quark is 
active, arranging combinations like (A c + A M ) 2 = A 2 which is numerically 



a We have studied this issue in much more detail in a simple model, with the same 
conclusion. 7 
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negligible since, unlike in the dimension-six case of Eq. (1), it cannot be 
compensated by a large m t factor. 

At n — m c , however, the charm quark is integrated out and the above 
mechanism no longer applies. At this scale, therefore, AS = 2 dimension- 
eight operators do get generated. Another simplification occurs, however. 
Notice that we eventually want to match onto the 0(p 2 ) chiral Lagrangian 
in Eq. (4) . Therefore any AS = 2 operator whose matrix element between 
a K° and a K is of higher chiral order is of no interest. Furthermore, in the 
large- N c limit four quark operators factorize, i.e. they become a product of 
two color singlets. A dimension-eight operator can be arranged in only two 
ways: either as a product of two dimension-four operators, or as a product 
of a dimension-3 operator times a dimension-5 one. However, the first case 
yields contributions of higher chiral order. This is due to the fact that 

(#%)|s l7i /ZV l |0) ~ , (21) 

which can be seen by contracting with g^ v and using the equations of motion 
-in the chiral limit-. So, only the combination C(dim — 3) x 0(dim — 5) is 
possible. However, sG IJ ' Vr y u d(x) is the only dimension-five current connect- 
ing a kaon to the vacuum 17 . Consequently, this means that there is only 
one AS = 2 dimension-eight operator to consider. This operator is given 

by 

Q2 _ 

Ce ff = C8 (^) 9sSLG^^d L (x) SL-f^dLix) +h.c. , (22) 

where g s is the strong coupling constant. 

Integrating out the charm quark yields for the coefficient c§ the matching 
condition 16 

c 8 (m c ) = - - \ 2 U - H A C A„. (23) 

O 

Below the charm mass the Effective Theory is given simply by 

Q2 _ Q 

Ceff = C 8(M) gsSLG^JvdL SL^f^dL - ~J=K SlJ^Ul ULjfidL + h.C. . 

(24) 

Notice that the second operator is the only one producing AS = 1 tran- 
sitions in the large- N c limit. Away from this limit its Wilson coefficient 
is not unity as in Eq. (24) but gets corrected by 30% — 40%, as a naive 
estimate of typical 1/N C corrections would say. 

At scales ji < m c , while one can still consider QCD in the perturbative 
regime, the Wilson coefficient c 8 (^) runs due to the fact that the square of 
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the AS = 1 operator in Eq. (24), with the two up-quarks propagating in a 
loop, mixes into the direct AS = 2 operator. At lowest order in the strong 
coupling constant, which is enough for our purposes, one finds 16 for scales 
/" ~ Aqcd ~ 1 GeV, 

u 2 

cg(n) = c 8 (m c ) + A 2 log — , (25) 

m c 

where c$(m c ) is given by Eq. (23). 

At scales fj, ~ Mk <C Aqcd perturbation theory is no longer a valid 
approximation and the matching of the Lagrangian (24) to the chiral La- 
grangian (4) requires again the machinery of the Hadronic Approximation 
to large- N c we used in the previous section. Skipping the details of this 
calculation 16 , we find that the coupling constant A| =2 gets the following 
contribution, 

A| =2 |dim-8 = jy A C A„ 5% + X 2 U (0.95 GeV) 2 j (1 ± 0.3) , (26) 

where S K is the parameter defined in Eq. (18), and the error is an estimate 
of the size of a typical 1/N C correction. Notice that all renormalization 
scale and scheme dependence has canceled out, as it should. This is one of 
the advantages of using the HA framework. 

The first contribution in (26) bears the mark of charm in A c and stems 
from the matching condition for the Wilson coefficient c 8 (m c ) when charm 
gets integrated out. Although it has an imaginary part that contributes to 
£k, its size is governed by the dynamical scale 5k ~ 350 MeV, which is 
small enough relative to m c to yield a very small correction to sk- 

The second contribution is the result of the running below the charm 
mass and subsequent matching onto the chiral Lagrangian (4); this is where 
resonances get integrated out. This is why only X u appears and, also, why 
the energy scale, which is a combination of resonance masses and couplings, 
is essentially given by the mass gap Aqcd ~ 1 GeV. Since X u is purely real 
it cannot contribute to Sk- The energy scale in (26) is comparable to m c , 
which explains why the Kl — K$ mass difference gets a sizeable <~ 10 — 20% 
correction 16 . 

3. e'/e and the A7 — 1/2 rule: "penguinology". 

At scales fi < m c the Standard Model gives rise to 10 four-quark operators 
capable of producing AS = 1 transitions 1 . All these operators mix as the 
renormalization scale evolves. 
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On the other hand, at the scale of the Kaon, the Effective Field Theory 
for AS = 1 transitions is described by the chiral Lagrangian 5 



Cfff 1 = --| V ud V* s {gsCs + 927 C 27 + e 2 g ew Tv (uX^U^Qr) } , (27) 



with = -iF 2 U{x)^D^U{x) , \ ( l 2) = 6 i3 6 j2 and Q L , R = 

diag(2/3,-l/3,-l/3). 



In this section we shall concentrate on the contribution to e'/e from Qe 
which, together with the one from Qs 27 ' 2S , is the dominant one 29 . As we 
shall see, the analysis will tell us interesting things about the AI —1/2 
rule as well. 

We may now run the Effective Field Theory from M\y down to the 
charm mass as usual but, at this point, make some simplifications in the 
analysis not to have to deal with the full 10 x 10 operator mixing matrix. 
Our first simplification will be to stay within the leading-log approximation. 
Furthermore we shall go to first subleading order in the 1/N C expansion but 
keeping only those terms which are enhanced by an extra factor of np 1 the 
number of flavors. In this case, the operator Q 6 only mixes with Q 4} where 



q—u,d.s q—u,d.s 

(28) 



In the previous expression the subscript MS is a reminder that these inte- 
grals are UV divergent and have to be regularized and renormalized using 
the same scheme as for the Wilson coefficients C4 j6 . This is also true for 
(ipip) and L5. As far as the A^ c counting goes, all the unfactorized contribu- 
tions, which arc of 0{np /N c ), are contained in the terms proportional to 
the functions Wdgrr and Wllrr- The terms proportional to L5 and unity 



where 




i=l,2,3 



and the sum over color indices within brackets is understood. 
The matching condition then reads 18 
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correspond to the factorized contribution from Qq and Q4 -respectively- 
and, formally, are of 0(N°). 

The functions Wdgrr and Wllrr are defined through the connected 
four-point Green's functions 

Wdgrm(<i) = A d'z (T{D Sq (x)G q MRUv)^s{z)}) , 

>«£f (?) = ^Jd A * dS d*z e*^T{I%(x)L&(0)R% u (y)R£.(z)}) , (30) 

after integration over the solid angle in g-momentum space, as in Eq. (8). 
In these expressions D Sq = s L q R , Gq d = q R dL, L% = s L ^q L , R^ u = 
d R j a uji, etc. It is not a coincidence that the pair of fermion bilinears 
which make up the operators Q6,4 in Eq. (28) also appear in these Green's 
functions, although they are located at different space-time points. It is 
the integral over Q in Eq. (29) which puts these two points back on top 
of each other. The matching condition (29) imposes that the same "mass 
term" r^J"™ is obtained when computed from the covariant derivatives in 
the Chiral Lagrangian (27) -which yields directly g$- as when computed 
from the four-quark Effective Lagrangian, which requires the insertion of 
the combination c 6 Q & + C4Q4 in the form shown in Eq. (29). 

In order to calculate the unfactorized contribution one should now con- 
struct the Hadronic Approximation. As in the previous case of Bk, an 
explicit cancelation of the rcnormalization scale dependence is achieved. 
This is good news since the scale dependence of the factorized contribu- 
tion is very large: it changes by a factor ~ 2 if the renormalization scale 
is varied in the range M p < \i < 1 GcV. Since this dependence has to be 
canceled by the unfactorized contributions, it is not unthinkable that these 
contributions be large. As discussed in the previous section on more general 
grounds, this could even be expected. 

And, indeed, it was found 18 that the unfactorized contribution from 
the operator Qq to the coupling g% in the matching condition (29) was a 
factor <~ 3 larger than the factorized contribution at a scale \i = 0.8 GeV. 
Although the effect was somewhat smaller for Q4, an enhancement was 
found there as well. Both large contributions come from what in the jargon 
is called an "eye" diagram (i.e. the contraction of the "dummy" quark q in 
the sums in Eq. (28)). On the basis of a Nambu-Jona-Lasinio model, this 
enhancement has also been found by Bijnens and Prades 20 . 

Decomposing the Wilson coefficients as 1 

c,(m)= Zi(fi) + Tyi(fi) , t= - * s td , (31) 
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the imaginary part of the coupling constant gs in Eq. (29) becomes 



where the error is the result of varying the quark condensate (which is 
the source of the biggest uncertainty), in the range (ip^p} 1 ' 3 (2 GeV) = 
(0.240 - 0.260) GeV. As to Imr, its current value 23 is -6.0(5) x 10~ 4 . 

We can also estimate Regs and 327- This is a very tough test for any 
calculational framework, as it is in these parameters that the AI = 1/2 
rule is rooted; a rule which still to this date defies detailed understanding. 
However, at the scale [i = m c , the situation is particularly simple since only 
the operators 



Q2 = ^{sLl^u L ){u Ll ^d L ) , Q x = 4(s L7 Ai d L )(u L 7 At w L ) (33) 



have nonvanishing real Wilson coefficients. Although at scales /x -C ra c pen- 
guin operators will again come into play, it is not crazy to neglect this effect 
and stay "as if fj, — m c " in a first approximation. In this case the contribu- 
tion to Regs from Qi j2 can be estimated because the non-eye diagrams are 
related to those appearing in the matching condition for gs=2 in the case of 
Bk, 25 while the eye diagram of Qi can be estimated from the eye diagram 
of Q4 in the matching condition (29), just setting nj? = 1. In numbers this 
leads to 18 Regs = 2.1 ± 0.8, to be compared to the experimental result 21 ' 24 
Regs = 3.6 ± 0.1, after subtraction of chiral corrections. In spite of the 
large errors involved I find this result quite encouraging, mainly because 
the strict large- N c result (factorization) would lead to Regs = 0.6 which is 
way too small. We begin to see the large unfactorized contributions which 
are indispensable for understanding the AI — 1/2 rule although, clearly, a 
more detailed analysis is needed before victory can be claimed. 

SU (3) allows to relate gs=2 in Bk to g 2 7 because their corresponding 
matching conditions to the four-quark Lagrangian are SU(3) rotations of 
each other 26 . This leads to 327 = gs=2 = 0.27 ± 0.12, to be compared to 
.927 = 0.30 ± 0.01 as extracted from experiment 21 ' 24 . 

Using all these results, together with 27 ' 28 



and, including final state interactions 30 and isospin breaking effects 21 , one 
obtains 18 e'/e = (5±3) x 10~ 3 to be compared to the experimental number 22 
(1.66±0.16) x 10~ 3 . This is not a bad result considering our approximations. 



lings — (3 ± 1) X Imr , 



(32) 



Im (e 2 g ew ) = (1.6 ± 0.4) • 10" 6 GeV 6 x Imr 



(34) 
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4. Application to lattice QCD: quenching penguins. 

Current calculations of e'/e done on the lattice require getting rid of the 
fermion determinant in the path integral in order to be numerically ef- 
ficient. This is accomplished by introducing some ghosts quarks which, 
although spin 1/2 particles, commute. This lattice technique, known as 
"quenching", has some dramatic consequences. In particular the flavor 
symmetry group is changed from the usual SU(3)l x SU(3)r to a graded 
SU(3\3)l x SU(3\3)b 31 ■ Furthermore, and even more importantly, there 
is no reason for the quenched theory to have the same weak low-energy 
constants as the true theory. The fact is that the current result for e'/e 
on the lattice is 3 times smaller than the experimental result, and with the 
opposite sign 32 . 

A particularly clear example of this is the transformation undergone by 
the strong penguin operator Qq in the quenched theory 33 . After quenching, 
the operator Q 6 in Eq. (28) is no longer a singlet under SU(3)r, or rather 
under SU(3\3)r- Instead, it can be decomposed as 

Q e = 2(32 7 X)(9fl7"«S + f R l^n) + 2(^ 7 X)(^7^ - feVSS) 

" v ' V v ' 

hQr q? ns 

(35) 

where the tilde on top refers to the ghost quarks. The flavor prop- 
erties are such that Q^ s is a singlet under the quenched flavor group 
SU (3\3)r whereas Qg NS is not. This changes their chiral representation 
accordingly 33 : 

Q® s => - ag' 15 STr(A£ M £^) + ag' 1} STr (2B A(UM + MU*)) 

Q% NS => f 2 a^ s STr (aUNU*) ; JV = \ diag (1, 1, 1, -1, -1, -1) (36) 

with qAq — sd, M the quark-mass matrix, and STr is the so-called super- 
trace. TV exhibits the non-singlet structure of Q^ NS . Comparing Eqs. (27) 
and (36), one sees that the couplings ag'^ have a counterpart in the true 
theory (the weak mass term of ag' 1 ^ is not written in Eq. (27)), but a^ s 
is a total quenching artifact. 

As we did in previous sections, it is now straightforward to apply large 
N c to determine these coupling constants a^' 1 ^, a^ s . One conclusion fol- 
lows immediately: because of the presence of the degenerate ghosts quarks, 
the sum over flavor in the internal propagator of the eye diagrams exactly 
vanishes. But this contribution was actually the dominant one in the QCD 
case!. Thus one obtains 34 : 
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(S 1) 

• the contribution from Q 6 to a q { ' is much smaller than its coun- 
terpart in the true theory, g$; and, 

a« s 1 

• = Terr; - 60 > 

where the quenched result 35 L 5 ~ 10~ 3 has been used. A recent lattice 
analysis 36 has obtained 

-jhr = 57(13) , (37) 

%i 

in nice agreement with our prediction (not postdiction!) above. 



5. Conclusions. 

Nature has an approximate chiral symmetry because the u, d, s quark 
masses are small. On a lattice, Nature is approached from the side of 
heavy quark masses by using chiral symmetry to guide the extrapolation. 
However, chiral symmetry is a property which is very difficult to achieve on 
the lattice, taking long hours of calculations with sophisticated codes and 
expensive computers and, regretfully, this becomes an important source of 
error. 

Large- iV c QCD offers the possibility to approach the problem from the 
other end. In the continuum chiral symmetry is much easier to achieve: 
it only takes the time needed to write a chiral Lagrangian. Furthermore, 
analytic calculations yield an understanding of the problem which allows 
building physical intuition. However, the flip side is that chiral symmetry 
is only exact when the quark masses are zero. To get to the real world, 
one must work one's way up to realistic quark masses by computing chiral 
corrections. 

I hope that the work presented here shows, among other things, how 
the continuum large- N c expansion may complement the lattice approach 
towards understanding kaon weak interactions. 
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